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We consider states with a charge- or spin-density wave order parameter which is odd in frequency, 
so that the order parameter vanishes at zero frequency and there is a conventional Fermi surface. 
Such states break translational symmetry and, therefore, are not conventional Fermi liquids. In the 
odd-frequency spin-density wave case, there are Goldstone bosons and the low-energy spectrum is 
manifestly different from that of a Fermi liquid. We discuss a simple model which gives rise to such 
ordered states. The frequency-dependence of the gap leads to an unusual temperature dependence 
for various thermodynamic and transport properties, notably the resistivity. 
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I. INTRODUCTION 



The Fermi liquid is a metallic state which is adia- 
batically continuable to a non-interacting electron gas. 
Its low-energy excitations are in one-to-one correspon- 
dence with those of the non-interacting system. Many 
attempts to uncover non-Fermi-liquid metallic behavior 
in strongly-correlated electron systems have focused on 
enhanced scattering mechanisms which might lead to 
anomalous behavior in the electron Green function. In 
this paper, we attempt a different tack and focus on 
metallic states with a conventional Fermi surface which 
can be distinguished from Fermi liquids by an order pa- 
rameter. In order to preserve the Fermi surface, the order 
parameter is taken to be odd in frequency |^ . When the 
order parameter breaks a continuous symmetry, the low- 
energy spectrum must also include Goldstone modes. As 
a result, odd- frequency ordering is a mechanism by which 
the set of low-energy excitations can be expanded from 
those of a Fermi liquid (so that they are no longer in one- 
to-one correspondence with those of a free Fermi gas), 
leading to a manifestly non-Fermi liquid state. 

A number of authors have considered superconducting 
states with order parameters {ca{k, co) cp{—k, —uj)) which 
are odd in frequency The energetic advantage of 

such a state is that it would enable the electrons to avoid 
instantaneous Coulomb repulsion while still benefitting 
from pairing. Some exactly soluble onc-dimensional mod- 
els, such as the one-dimensional Kondo lattice, have a 
tendency towards such ordering |8|j^. However, there 
are claims that the simplest models of odd-frequency su- 
perconducting states suffer from pathologies which make 
them unstable |[7|,p^. 

In this paper, we consider the analogous states in 
the particle-hole channels. We find that simple mod- 
els of odd-frequency density wave states do not suffer 
from any pathologies (and suspect that this holds for 



the corresponding superconducting states, as well). Odd- 
frequency states can exhibit a number of interesting non- 
Fermi liquid properties including the Goldstone modes 
mentioned above, a range of states above the Fermi sur- 
face with finite lifetimes even in the limit of vanishing 
frequency and temperature e,T —t 0, and a non-mean- 
field-like temperature-dependent order parameter. 

However, the question of their detection is non-trivial 
and cuts to the heart of attempts to experimentally dis- 
tinguish non- Fermi liquids from Fermi liquids. While an 
odd-frequency superconducting state is, first and fore- 
most, a superconducting state, which would be identi- 
fied by its vanishing resistivity, Meissner effect, etc., an 
odd-frequency density wave state can masquerade as a 
Fermi liquid since the order parameter vanishes at zero 
frequency. There will be signatures in thermodynamic 
and transport measurements, but they can easily be mis- 
takenly attributed to other effects, as we discuss below. 



II. ORDER PARAMETERS AND SYMMETRIES 

An odd-frequency charge-density wave state is defined 
by the anomalous correlation function 

(c"t (fc, e„) (fc + Q, e„)) = F {k, e„) , (1) 

where F (fc, e„) is an odd function of frequency, 

Fik,-en)^-Fik,en). (2) 

To find an equal-time correlation function which serves 
as an order-parameter, we Fourier transform ([^): 



{Tr (c"t {k,T)c^ik + Q,0)))^Fik, t) 



(3) 



where F (fc, r) is imaginary and odd in r. Hence, we can 
use the time-derivative as an order parameter: 

{Tr (9.c"t {k, t) C„ (fc + Q, 0)))^^„ = drF {k, 0) . (4) 
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The state defined by these order parameters, (|l|) or (|^), 
breaks translational symmetry. Time-reversal symmetry 
is not broken. This is most easily seen by considering (||). 
Taking the complex conjugate of both sides of (|^) gives 

{Tr (c"t (fc + Q, 0) c„ (A;, r)) ) - (f (fc, r)^ 

^-F{k,T). (5) 

Meanwhile, (|^) is transformed under time-reversal, T, 
into: 

{Tr{T (c"t (fc,r) c„ (fc + Q,0)))) 

= (r,(c"t (fc + Q,0)c„ (fc,-T))) 

= -F(fc,-T) 

= (T,(c"t(fc,r)c„(fc + g,0))), (6) 

and, hence, the order parameter does not break time- 
reversal symmetry. In going from the first equality to 
the second, we have used (||). 

An odd-frequency spin-density wave is defined by 

(c"^ (fc, e„) cp {k + Q, £„)) =ri-d$F (fc, e„) , (7) 

where F (k, e„) is again an odd function of frequency and 
n is the direction chosen spontaneously by the ordered 
state. This state breaks translational symmetry and spin- 
rotational symmetry, which is broken to the U(l) sub- 
group of rotations about n. Again, time-reversal is pre- 
served. Such a state will exhibit a non-zero expectations 
value and anomalous correlations of the spin nematic or- 
der parameter, SiSj — SijS'^/3, such as those discussed 
for spin-only models in uM- 

III. MODEL INTERACTION 

We now consider a simple two-dimensional model 
which admits an odd-frequency charge-density-wave 
state at the mean-field level. The model contains a non- 
singular four-fermion interaction which can be generated 
by the exchange of phonons or some gapped electronic 
collective mode. For simplicity, we focus on the charge- 
density-wave case; the spin-density-wave is analogous. 

We consider an effective action which consists of a ki- 
netic term 

So^ fdr f 4^ c"t (fc, r) (5, - (e(fc) - a^)) c„ (fc, r) 



(2^)2 



(8) 



and an interaction term 
1 



'S'int — 



dr 



k,k' 



e\k + Q)drCo,{k) 



Vkk' 



{k')drCp{k' + Q) 



(9) 



To avoid clutter, the J^, j,, is used as shorthand for the 
integrals over k,k' . The subscript c indicates that the 
terms in brackets are actually defined with a frequency 
cutoff, ilr- 



C"\k + Q)drCjk) 



{k,en). (10) 

\e„\<^a 

For simplicity, we take Vkk' independent of k,k' , Vkk' — 
A. For simplicity, we also take Q = (tt, tt) and e(fc) = 
— 2t(cos kx + cos ky), corresponding to commensurate or- 
der for a system of electrons on a square lattice with 
nearest-neighbor hopping. The generalization to incom- 
mensurate order and other band structures is straight- 
forward. 

The interaction term S'int is long-ranged in precisely 
the same way as the BCS reduced interaction. A more 
realistic short-ranged interaction would be of the form 



Sint — 



k^k' ,q 
X 



C'^^k + q)drCa{k) 
C^\k')drCp{k' + q) 



v 



(11) 



which includes (^) as one term in the sum over q. Such 
an interaction could arise from the diagram of figure ^ if 
the collective mode has a propagator of the form 



nl + v-'ik-k'Y + ien-en'f 



(12) 



If V is small, then we can expand the collective mode 
propagator to obtain Vkk' = A, for |e„| , |e„'| < Q,c- Other 
terms will also be generated which could drive the for- 
mation of even-frequency order, but they appear to be 
weaker. 

k+q,& + (i) k'+q,z'+(i) 




k\ e' k, e 

FIG. 1. A diagram which can lead to an interaction fa- 
vorable for odd-frequency density-wave ordering. The dashed 
line represents a collective mode which mediates the interac- 
tion. 
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We assume an order parameter of the form: 
A f (Pk 



(c"t(fc + Q)0,c„(fc))^ 



(13) 



According to our previous observations, a is real. Then 
the mean-field action takes the form 

-5*0 = / c"'^(fc, e„) (ie„ - (e(/c) - /i)) Ca(fc, e„) 



(2^ 



Q!e„c°'l'(fc,e„)ca(fc + (5,e„). (14) 



The equation of motion following from the mean-field ac- 
tion is: 

(ie„ - e(fc) + n)ca{k,en) - aenCa{k + (9,e„) = 0. (15) 

We multiply the equation of motion by cjj(fc, e„) and take 
the imaginary time-ordered expectation value. We see 
that the ordinary and anomalous Green functions satisfy 
the equation: 

(ie„ - e(k) + G{k, e„) - ae„i^(fc, e„) = 1. (16) 

The right-hand-side results from the time-derivative (i.e., 
«e„) acting on the time-ordering symbol p2) . If we make 
the replacement k ^ k + Q va. (|l^) , then we can derive a 
second equation in the same way, 

(ie„ -f e{k) + fi) F{k, e„) - ae„G(fc, e„) = 0, (17) 

and the Green function and anomalous Green function 
for |e„| < n,c take the form: 



G(fc,e„) 
F(fc,e„) 



ien + fi + e{k) 



(ie„ + Ai) - {e{k)) - a'^el 

acn 

(ie„ + fif - {t(k)f - a2e2 



(18a) 
(18b) 



One might naively think that equation (y_7|) could be ob- 
tained by inspection from (|6|) by replacing A(e„) by 
A*(— e„) as one usually does in the case of an even- 
frequency gap. In this case, this would amount to the 
replacement of ae„ by — ae„, as was done in Ref. ||7|,|lO|. 
However, in the case of an odd-frequency gap, this sim- 
ple substitution only works for reaZ frequencies. Since the 
gap is linear in frequency and the squared modulus of the 
gap appears in the Green functions, the Green functions 
are no longer analytic in the frequency. As a result, the 
analytic continuation from real to Matsubara frequencies 
is subtle. The safe route is to derive both (16) and ( |l^ ) 
directly from the mean-field action, as we have done. The 
naive, incorrect form of (^|) would lead to a negative su- 
perfluid density in the case of odd-gap superconductors. 
On the other hand, the correct analogue of ( p^ ) for a 
superconducting action induced by disorder was used in 
Ref. @. 



For |e„| > fic, vanishes and G returns to its normal 
state form. In principle, we should also allow for a quasi- 
particle renormalization Z resulting from the interaction, 
but this does not qualitatively modify our results, so we 
drop this correction for simplicity. 

From these Green functions, we see that, at /i = 0, 
the odd-frequency charge-density-wave order parameter 
modifies the quasiparticle spectrum to: 



(19) 



i.e., it renormalizes the effective mass. For /i ^ 0, the ef- 
fect is more complicated. As a result of the odd-frequency 
charge-density-wave, e(fc) — /i is replaced with 



Eik) . -^-^^ 



a2) e2 (fc) 



1 



(20) 



From (|2y), we see that the Fermi surface is unmoved, i.e., 
E{k) = when e(fc) — ji, as we expect, since the order 
parameter vanishes at zero frequency. 

Furthermore, there is a range of k values above the 
Fermi surface, |e(A;)| < [a/i/Vl + Q;^ | , where E(k) has an 
imaginary part, so that quasiparticles in this region have 
a finite lifetime. However, these states have zero occu- 
pation number, as the corresponding poles in the Green 
functions turn out to be lying outside of the integration 
contour. This gives rise to a possibility of the ground 
state in which some of the quasiparticles are in levels 
which are in disconnected from the rest of the Fermi sea. 



IV. GAP EQUATION 

We must now impose a self-consistency condition on 
F (fc, e„), which is the gap equation. We will also impose 
a condition on the particle number, thereby implicitly 
fixing the chemical potential. These conditions read: 



"c „ Jkiie 



^Jk{ien+fi) - (e(fc)) -a^el 
yV^' f ^e„ + /i + e{k) 

„ Jk{ien+ fJ.f - {e{k)f ~ a^el 



1, (21a) 
n. (21b) 



For simplicity, we consider the case of half-filling, n ~ 1. 
We have repeated our calculations at non-zero doping 
and found similar results. 



The prime on the Matsubara summations in (21a), 
( [21 hp indicate that they are done with a = for |e„| > flc 
and a only for |e„| < ilc- A more realistic model 
replaces the interaction with one that has a "smooth" 
cutoff Sriitn) Sri{tn') and a by as,,(e„), with s,,(e) = 1 



for e ^ and Sj,(e) 
•q between the limit rj 



for e 3> ric- We can vary 
> oo, which corresponds to a 
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sharp cutoff, and 77 ^ which corresponds to the ab- 
sence of a cutoff. For computational simpUcity, we take 
Sr,{e) = riF(|e| - ^Ic, P = v)- 

We now discuss the analysis of (21a), (21b). Let us 
first consider the case of a sharp cutoff. The left-hand 
side of the gap equation vanishes if the temperature is 
above 



(22) 



jimax ^j^g highest possible transition temperature for 
this system. Just below this temperature there is only 
one pair of terms in the Matsubara sum which is allowed 
by the cutoff fie. As the temperature is decreased, more 
Matsubara frequencies begin to contribute, resulting in 
minor steps in the phase diagram. 

For large A, Tc — T^^'^^- Decreasing A, we enter a 
regime, Ac2 < A < Ad, in which the system is in the odd- 
frequency density-wave phase for an intermediate range 
of temperatures Tc2 < T < Td- Ad is the location of 
the quantum phase transition at which the odd-frequency 
density- wave orde r firs t appears at zero temperature. It 
is obtained from ( ^la| ) by setting /i = and a = and 
converting the Matsubara sum into an integral: 



(2^ 



fic — e {k) arctan 



e(fc) 



Ad 



(23) 



The integrand is approximately for small e (fc) and 
n3/3e2 (fc) for |e(fc)| > n^. 

As A is further decreased, Tc2 increases and finally 
reaches T™'^'^ at Ac2. To find Ac2, we again set ^ = 
and a = 0, but now we retain only the pair of terms 
in the Matsubara sum corresponding to the frequencies 



/ 



(fk 



{2nf + e2 (fc) 



1 

Ac2 



(24) 



Comparing expressions (|23| ) and (|2j), we find that the 
latter is larger by a factor of 3 for large e(fc); hence, 
Ac2 < Ad. For A > Ac2, a jumps discontinuously at 
j.max^ Finally, for < A < Ac2, there are no odd- 
frequency charge-density-wave solutions. As a result, the 
phase diagram has the shape shown in Fig. || (dashed 
line) with re-entrant transitions for Ac2 < A < Ad. 

Let us now consider how this picture is modified when 
we make the cutoff smooth, as it must be in a physical 
system. The sharpness of the steps which separate the 
re-entrant transitions depends on the details of the high- 
frequency cutoff; they disappear in the limit that the cut- 
off is very smooth. When the cutoff is relatively sharp, 
more re-entrant transitions are possible, and there will be 
several temperature regions in which an odd-frequency 
charge-density wave occurs. However, for a smooth cut- 
off, there is typically only one such region. The transition 



at T^'^^ is replaced with a smooth curve Td(A), at which 
a second-order transition takes place. The corresponding 
phase boundary is depicted by the solid line in Fig. |[ 




Normal phase 



0.00 0.02 0.04 0.06 0.08 0.10 

T T^" 
FIG. 2. Pliase diagram at half-filling for "smooth" and 
"sharp" cutoffs at = 0.3. The energy scale is set by the 
width of the band W = 2. 



Below rd(A), a increases as shown in Fig. |[ Note 
that for a smooth cutoff all transitions are of second or- 
der, even though the rise of a at Td is very steep for large 
values of A and may give a false impression of a first-order 
transition. It is also noteworthy that even for large A the 
order parameter a rapidly attains its maximum as T is 
decreased below Tc and then decreases as T — > to some 
non-zero asymptotic value. This has a significant impact 
on experimentally measurable parameters, as we describe 
in the following section. 




FIG. 3. Order parameter for different values of the inter- 
action strength. A, with a "smooth" cutoff at = 0.3. 

The unusual temperature dependence of the order pa- 
rameter, which is reflected in the re-entrant phase dia- 
gram and (as we will see in the next section) the con- 
densation energy, is a consequence of the frequency- 
dependence of the gap. In an ordinary, frequency- 
independent (or weakly-dependent) ordered state, the 
condensation energy at weak-coupling comes primarily 
from states near the Fermi energy. At high tempera- 
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tures, these states are thermally excited, so there is little 
condensation energy to be gained, and the order param- 
eter decreases as the temperature is increased. In the 
case of odd-frequency ordered states, there is very lit- 
tle condensation energy to be gained from the particles 
near the Fermi surface because their energy is low (and, 
hence, they interact weakly with the order parameter). 
As a result, the order parameter and condensation energy 
decrease as the temperature is decreased. 

In fact, there is a second solution to the gap equation in 
which it is favorable to have a disconnected Fermi sea in 
which some electrons are excited to a strip in momentum 
space above the Fermi surface which is diconnected from 
the rest of the Fermi sea, which is centered about k = 0. 
When this occurs, the occupation number does not in- 
crease monotonically with e{k). However, this solution is 
higher in energy, so it does not occur. 



V. EXPERIMENTAL SIGNATURES 

At Tc, there will be the usual thermodynamic signa- 
tures of a second order phase transition. The condensa- 
tion energy associated with this transition is the differ- 
ence between the free energy of the odd-frequency charge- 
density-wave state and corresponding free energy of the 
normal state at the same temperature. At fi = and 
with a sharp cutoff, 



- in a complete reversal of the situation for a frequency- 
independent gap ]T^ . 

Again, for a smooth cutoff the Matsubara frequency 
sum becomes infinite and a should be replaced with 
s^(e„)a. For smooth cutoff and finite temperature, the 
condensation energy must be evaluated numerically. The 
dependence of the condensation energy AE on tempera- 
ture is shown in Fig. ^. Unlike in even- frequency phases, 
where \AE\ slowly increases as the system cools down 
and attains its maximum at zero temperature, in our 
model it rapidly reaches a maximum and decreases to a 
constant asymptotic value as T ^ 0. As may be seen 
in figure ^, the condensation energy is of order of NpT^i 
at the maximum, which is comparable to the maximum 
condensation energy attained in even-frequency phases. 



< 




AE{T) = 



2T 



E 



In 



(1 



(25) 



where — {Qc/t^T — 1) /2. This equation is obtained 
by using a as a Hubbard-Stratonovich field to decouple 
(^, resulting in the first term in (|25|). The electronic 
action is then the mean- field action (nj) , so that the par- 
tition function may be evaluated to give the second term 
in (25). At zero temperature, the integrals may be eval- 
uated analytically: 

AE(0) ^-ra^ +- e(fc) arctan — ^ 



e{k) 



T^JkVT+a 



arctan 



2 ^ e{k) 

Og + £^(fc) 



In 



(26) 



The last term is overwhelmingly negative, as may be seen 
in various limits (e.g., a ^ 1 or ilc oo). Note that the 
energetic gain comes not from the terms in the frequency 
sum ( p5| ) with small Matsubara frequency, which actu- 
ally increase energy, but from the terms near the cutoff 



FIG. 4. Condensation energy as a function of temperature 
in model (0) with a "smooth" cutoff at fl — 0.3. 



The DC conductivity can be computed from the Kubo 
formula. We assume a model in which impurities give 
rise to a lifetime, r. The resulting conductivity is given 
by 



dejk) 
dkx 



AT 



cosh 



E{k) 
2T 



(27) 



where E{k) is given in (|20|). For small /i, the modifica- 
tion of the quasiparticle spectrum in the odd-frequency 
charge-density-wave phase reduces to the rescaling of the 
electron mass near the Fermi surface, so that the new 
effective mass is m* = m^/l + o? . This mass enhance- 
ment leads to a noticeable increase of the resistivity, as 
shown in Fig. |^. Of course, outside of the the region 
7c2(A) < T < Tci(A) the resistivity is that of the normal 
phase. 
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FIG. 5. Resistivity as a function of temperature in the 
model of Eq. (^) with a "smooth" cutoff a,t Q — 0.3 and 
interaction strength A = 4. 



VI. DISCUSSION 

In this work we have considered a model with an odd- 
frequency charge-density-wave solution. The transition 
to this state is signaled by a second-order phase transi- 
tion with a jump in the specific heat. For strong interac- 
tions, the model is in an odd- frequency charge-density- 
wave phase for all temperatures T < Tc- For moderately 
weak interactions, the model is in such a phase for an 
intermediate temperature regime Td < T < Tc2- There 
is a quantum phase transition at A = Aci; for weaker 
interactions, the system is not ordered at T = 0. 

A similar model admits an odd-frequency spin-density- 
wave ground state. Such a state will have, in addition to 
its Fermi-liquid-like quasiparticles, Goldstone boson exci- 
tations. As a result SiSj — SijS'^/3, has a non-zero expec- 
tation value, and its correlation functions have Goldstone 
poles. 

The method that we used to derive the Green functions 
of the odd-frequency density wave can be applied to the 
analysis of odd-frequency superconductivity as well, with 
some modification. The corresponding mean field theory 
should describe a stable state with positive superfluid 
density. 

Odd-frequency charge-density wave order results in 
mass enhancement. This affects transport properties; 
in the density-wave state the resistivity is considerably 
larger than in the normal state. The effect is largest 
at intermediate temperatures. At low temperatures, the 
system is either in the normal state (for A < Ad) or 
in an ordered state (for A > Ad) with some asymp- 
totic value of the order parameter. This type of non- 
monotonic resistivity curve has been observed in a num- 
ber of strongly-correlated electron systems. In layered 
materials such as the cuprates and ruthenates, it has 
been observed in c-axis transport In 2DEGs, this 

type of behavior has been observed in the vicinity of a 
putative metal-insulator transition . It would be pre- 



mature to suggest that odd-frequency order is developing 
in any of these experiments, but it is noteworthy that it 
does provide a natural explanation of otherwise puzzling 
behavior. 

Odd-frequency density wave ordering is also mani- 
fested in thermodynamics. Unlike in even-frequency 
states, where the condensation energy \AE\ is small near 
the phase transition and reaches a maximum at zero tem- 
perature, in odd- frequency states the maximum of |Aii^| 
is located near the upper critical temperature. Conse- 
quently, there is strong variation of all thermodynamic 
quantities with temperature just below Tc- Again, at 
lower temperatures these phenomena disappear. 

The model that we have introduced is certainly simpli- 
fied, as it ignores the possible proximity of other phases. 
However, we believe that odd-frequency density-wave or- 
der can result when a tendency towards ordinary even- 
frequency density-wave order is frustrated by some com- 
peting interaction. The resolution of the competition will 
depend on the scales at which the various interactions 
act - essentially ilc in our model. In particular, one can 
imagine a scenario in which an ordinary even-frequency 
density wave is favorable at higher temperatures, but be- 
low a certain temperature, the system undergoes a tran- 
sition into an odd-frequency state. 
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